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Abstract 

We investigate the possibilities of searching for non-standard CP violation in pp —>■ 
ttX at multiple TeV collision energies. A general kinematic analysis of the underlying 
partonic production processes gg — > tt and qq tt in terms of their density matrices 
is given. We evaluate the CP-violating parts of these matrices in two-Higgs doublet 
extensions of the standard model (SM) and give results for CP asymmetries at the 
parton level. We show that these asymmetries can be traced by measuring suitable 
observables constructed from energies and momenta of the decay products of t and 
t. We find CP-violating effects to be of the order of 10~^ and show that possible 
contaminations induced by SM interactions are savely below the expected signals. 
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1 Introduction 



A high energy and high luminosity proton-proton colhder, such as the planned Large 
Hadron Collider (LHC) at CERN, would be capable of producing millions of top and 
antitop quarks. This would offer the unique possiblity to explore in detail the physics 
of these quarks - which have not been discovered yet, but for whose existence there is 
indirect evidence |T|]. Specifically, since the top is known to be heavy, rrit > 113 GeV 
[0], precision studies based on large samples of ti events may serve as a probe, through 
sizeable top- Yukawa couplings, to the electroweak symmetry breaking sector (Higgs 
sector for short). This sector may have a richer structure than the one conceived in the 
standard model (SM) — as is the case in many of its extensions. As a consequence a 
number of new phenomena may exist. A particularly intriguing one is a new "source" 
of CP violation provided by the Higgs sectoiQ which is unrelated to the Kobayashi- 
Maskawa phase This is possible already in the two-Higgs doublet extensions of 
the SM 0-0: here neutral Higgs boson exchange leads to CP- violating effects in 
fermionic amplitudes, and these effects would show up most pronouncedly in reactions 
involving top quarks The subject of this paper is to investigate in detailQthe 



manifestation and the magnitude of neutral Higgs particle CP violation in pp ttX. 
(For other studies on CP violation in top quark production and decay see |p!2|-p5|.) 



The outline of our paper is as follows: In section 2 we give a general kinematic 
analysis of the reactions gg —>■ tt and qq — *■ tt which are the leading partonic processes 
in pp —y ttX. We study these reactions in terms of their production density matrices 
and describe the properties of these matrices under various symmetry transformations 
including CP transformations. In section 3 we evaluate the CP-violating parts of the 
density matrices in a specific model, namely the two-Higgs doublet extensions of 
the SM with CP-nonconserving neutral Higgs boson exchange. Correlations which 
are sensitive to CP violation at the parton level are identified and results for their 
expectation values are presented. In section 4 we show that these CP asymmetries 
can be traced in pp — > ttX by looking at simple observables which involve energies 
and/or momenta of the decay products of t and t. Moreover, possible contaminations 
by CP-conserving interactions are discussed and shown to be much smaller than the 
expected signals. In an appendix we list the analytic results of our calculations. 



"'^This source was shown to be of interest in attempts to explain the cosmological baryon asym- 
metry g. 

short account of our work was given in 



2 Production density matrices for gg tt and qq 

ti 

Because a heavy top has an extremely short hfetime (r^ < lO^^^s if nit > 100 GeV), 
the polarization of and spin-spin correlations between t and t are not severely diluted 
by hadronization |^ . These are "good" observables in the sense that effects involving 



the spins of t and i can be treated perturbatively. Therefore we will discuss the reac- 
tion pp — *■ ttX, respectively the underlying partonic processes in terms of production 
density matrices. We will consider only unpolarized pp collisions. At LHC energies tt 
pairs are produced mainly by gluon gluon fusion. This reaction dominates over quark 
antiquark annihilation into ti. We will first discuss the (unnormalized) production 
density matrix for the reaction g{pi) +g{p2) ^(^i) +^(^2) in the gluon-gluon center 
of mass system. It is defined by 

Rla,, A/5.(P, k) = N^' E {tiki, t(A;2, Pi)\r\gip,), gip,))* 

colors, gluon spin 

{tih,a2),i{k2,P2)\r\gip,),gip2)) (2.1) 

where a, P are spin indices, p = Pi, k = ki and Ng = 256. We sum here over the 
gluon spins and colors since we are only interested in analysing the polarizations of 
the t and i and their spin-spin correlations. The matrix can be decomposed in 
the spin spaces of t and i as follows: 



= An^l + Bf^a' ® U + Bf-1 ® a' + Cf^a' ® . (2.2) 

The first (second) factor in the tensor products of the 2x2 unit matrix U and of the 
Pauli matrices a* refers to the t (t) spin space. 

Because of rotational invariance, the functions Bf^ and Cij can be further decom- 
posed: 



Cij = CgoSij + eiji{cgipe + Cg2ke + Cg^hf) 
+Cg4PiPj + Cg^kikj + Cgeipikj + pjki) 
+Cg7{pinj + Pjhi) + Cgsikfij + kjhi). (2.3) 

Here the hat denotes a unit vector and n = p x k. The structure functions A^, and 

Cgi depend only on s = (pi +^2)^ and on the cosine of the scattering angle, z = p ■ k. 

Next we discuss the properties of under various symmetry transformations. 
Since the initial gg state is Bose symmetric, R^ must satisfy 



2 



it:^(-p,k) = i?^(p,k). 



(2.4) 



The initial gg state, when averaged over colors and spins, is a CP eigenstate in its 
center of mass system. It is therefore possible to classify the individual terms in 
according to their CP transformation properties. If the interactions were CP- 
invariant, the matrix would have to satisfy 

RU,a,aM^) =<a„/3./5.(P,k). (2.5) 

In table 1 wc give a complete list of the transformation properties of the structure 
functions under P, CP, and exchange of the initial gluons (''Bose"). It is also in- 
structive to collect the properties of these functions under time reversal (T) and CPT 
transformations neglecting, just for this purpose, non-hermitean parts of the scatter- 
ing matrix. To give an example, table 1 is then to be read as follows: for a T-invariant 
interaction one has bf^{z) = —bf^{z), i.e., 6^3 = only at the Born level, whereas at 
higher orders absorptive parts render this function non-zero. 

Because of Bose symmetry, the structure functions A^, bf2, Cgo, Cg2, Cg^, Cg^ and 0^7 
are even functions of z, the other functions are odd in z. 

The contributions to R^ can be decomposed into a CP-even and a CP-odd part: 

= Rlen + R'CP- (2.6) 

As can be read off from table 1, the CP-even term R^-^^^ in general has the following 
structure: 

^?ven = ^^11 ® 11 + (6 + 6^^"^^, + 6^^"^n,) ((7^ ® 11 + 11 ® a') 
+(cgo% + CgiPiPj + Cg^kikj + CgQ^pikj +pjki) 
+Cg7{pinj + pjhi) CgsikiUj + kjhi))a^ ® aK (2.7) 

Nonzero Ifgi"^, c^y, c^g can be induced only by parity-violating interactions, 

C(,7, Cg^ need in addition absorptive parts in the scattering amplitude when the in- 
teractions are CPT-invariant. The structure function 6^™" can only get contributions 
from absorptive parts induced by parity-invariant interactions. 
The CP-odd term i?^p reads 

Rcp = {b^giPi + b^'gik + h^Pfii) (a^ ® 11 - 11 ® 

+^ijkicgiPi + Cg2ki + Cg3ni)a^ (8) a''. (2.8) 
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CP-violating interactions which are also parity-violating can give contributions to 
h^i , b'^^ , Cgi, Cg2- Nonzero fo^/', 6^2^ require in addition absorptive parts. C- and CP- 
violating interactions can induce non- vanishing structure functions b'^^, Cgs, where 
Cg3 7^ requires in addition absorptive parts. 

The above discussion of the transformation properties of the structure functions 
holds to all orders of pertubation theory. 

The production density matrix W for qq ti is defined in complete analogy to 
O as 

A/5.(P,k) =iV-i Yl m,a^),iih,(3,)\T\qip,),qip2)y 

colors, gq spins 

{t{h, a2),i{k2, P2)\r\qipi), q{p2)), (2.9) 

where A'^, = 36. The decomposition of in the spin spaces of t and i is exactly the 
same as for (equ. (|2.2|) , (|2.3| )) as is the splitting into CP-even and odd terms (equ. 
( |2.6| )- ( p.8| )). The transformation properties of the structure functions A'', . . . , Cgs of 
are the same as the respective ones for A^, . . . ,Cgs of given in table 1. Thus all 
conclusions derived from these transformation properties — except those from Bose 
symmetry, of course — are also valid for the structure functions of R*^. 

3 CP violation and density matrices in two-Higgs 
doublet models 

Up to now our discussion has been independent of any specific model. Suffice it to 
say that the Kobayashi-Maskawa mechanism of CP violation ^ induces only tiny 
effects in the flavor-diagonal reactions of sect. 2. In the following we will concentrate 
on CP-violating effects generated by two-Higgs doublet extensions of the SM with 
CP violation in the scalar potential . We briefly recall the features of these models 
relevant for us. CP violation in the scalar potential induces mixing of CP-even and 
-odd scalars, thus leading to three physical mass eigenstates \ipj) (j = 1, 2, 3) with no 
definite CP parity. That means, these bosons couple both to scalar and pseudoscalar 
fermionic currents. For the top quark these couplings are (in the notation of [^): 

3 

Cy = -{V2GfY^^ Y^{ajtmtit + djtmtiij5t)ipj, (3.1) 

where Gp is Fermi's constant, rrit is the top mass, 

ajt = c?2j/ sin /5, djt = —d^j cot /3, (3.2) 

tan/? = V2/V1 is the ratio of vacuum expectation values of the two doublets, and 
(i2j, d^j are the matrix elements of a 3 x 3 orthogonal matrix which describes the 
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mixing of the neutral states [Q. In the following we assume that the couplings and 
masses of V92,3 are such that their effect on all quantities discussed below is negligible. 
Then the measure of CP violation generated hj ip = ipi exchange in flavor-diagonal 
reactions like gg tt, qq ti is 

■jcp = ~<^0j = d2id3i cot P I sin /3, (3.3) 

where we have put a = au, d = du- So far, data from low energy phenomenology, 
in particular the experimental upper bounds on the electric dipole moments of the 
neutron and of the electron do not severely constrain this parameter: ^cp 



may be of order one. We note here that the couplings of the ip^ to quarks and leptons 
induce CP violation already at the Born level. The especially interesting case of a 



Higgs boson Lp decaying into tt was shown in to lead to CP- violating spin-spin 



correlations which may be as large as 0.5. (For other discussions of the CP properties 



of neutral Higgs bosons see [BS 



We will now discuss the structure of the matrices i?^p and -R^p in these models. 
The Higgs boson contributions to the processes gg ti and qq ti discussed in 
section 2 are shown, together with the leading SM diagrams, in figs. 1, 2. Since 
the CP-nonconserving neutral Higgs exchange is, in particular, parity-violating, the 
relations 

i?^j^Vp,-k) = -i??j?(p,k) (3.4) 

hold as long as Rcp results from interference of these Higgs exchange amplitudes 
with amplitudes from parity-invariant interactions. This forces bgf, Cgs and b^f, Cqs 
to be zero in these models. Furthermore, the virtual intermediate gluon produced 
by annihilation of unpolarized q and q cannot have a vector polarization. Thus the 
contributions of fig. 2 to are invariant with respect to the substitution p = Pi — > 
—p. Hence the structures of R^p and R^p are the same; the functions bgf, Cgi, b^f, 
Cqi of eqn. ( p.8|) are odd under z —z, whereas bg^, Cg2, b^^f, Cq2 are even functions 
of z. 

The explicit results for the matrices R^ and R'^ evaluated from the diagrams of fig. 
1, 2, respectively, are given in the appendix. The width of ip must be taken into 
account in the calculation of R^ ii ip > 2mt, since in that case the contribution from 
fig. Ih can become resonant. Because for (p > 2mt the width of (p is not very small 
as compared to its mass it is important to note that the narrow width approximation 
cannot be applied in this case. 

In view of the above discussion it is now very easy to identify the correlations at the 
parton level which trace the various CP-odd parts of the production density matrices. 
The expectation value of an observable O for the respective parton reactions is defined 

as 

j\dztr(R'0) , , , , 

{O), = , ^ — — (i = g, q) 3.5 



5 



Contributions of the functions h'^ig2i b^iq2 picked up by taking expectation values 
of 

k-(s+-s_)/,(z), (3.6) 

or 

p-(s+-s_)/o(z), (3.7) 

or linear combinations thereof, where s+, s_ are the spin operators of t and t, respec- 
tively, fe{z) is an even function of z and fo{z) is odd in z. One has for example 



(k ■ (s - — s_))o = T—^ — , (3. 



and likewise for (k ■ (s+ — s_))g. 

The result for the basic longitudinal polarization asymmetry (k ■ (s+ — s_))g is 
plotted in fig. 3 as a function of the parton CM energy for rrit = 150 GeV and 
two values of the Higgs boson mass: The dashed curve corresponds to = 100 
GeV and 'jcp = 1- For "^.^ order of 2mj or larger, the shape of the resulting 

graph depends, for fixed m^, on the strength of the Higgs couplings a, a and on 
the couplings to W~^W~, ZZ determining the width of the ip. (See eqn. ( [A.10| ) for 



details.) For the sohd line we have chosen = 350 GeV, \a\ = \a\ = •jcp = 1 
and = 47 GeV. The asymmetry (k ■ (s_ — s+)) corresponds to the asymmetry 
ANlr = [N{tLiL) - N{tRiR)]/{a\\ tt) studied in 0. We reproduce the numerical 
results of for AN^r if we neglect and use 7cp = 1/v^ which corresponds to 
the parameter Im(A^) = ^/2 used in |]10 . 



The functions Cgi^g2, Cgi_g2 generate nonzero expectation values of the triple prod- 
uct correlations 

k- (s+ X s_)/ie(z), (3.9) 

p- (s+ X s_)/i„(z), (3.10) 

and of their linear combinations, where hg and ho are even and odd functions of z, 
respectively. For example. 

In fig. 4 we plot this basic CP-odd and T-odd spin-spin correlation with the same 
choice of parameters as in fig. 3. It reaches values of up to about two percent. 

For completeness, we show in figs. 5 and 6 the expectation values (k • (s+ — s_))q 
and (k ■ (s+ x s_))g, respectively, again for the same choice of parameters as in fig. 
3. Here the CP asymmetries get smaller with growing Higgs masses. 
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4 CP observables for pp ttX 

The CP-violating spin-momentum correlations for t and t of the previous section must 
be traced in the final states into which t and i decay. In this section we discuss a 
few observables which allow to do this. The charged lepton from t Wb — ^ i~^i>ib is 
an efficient analyzer of the top spin [^. We will therefore conider only decay chains 



where at least one of the top quarks decays semileptonically. We shall use the SM 



decay density matrices as given in [16, 1 



Observables in pp tiX cannot be classified as being even or odd with respect to 
CP, because the initial state is not a CP eigenstate. However, they can be classified 
as being T-even or T-odd (i.e. even or odd under refiection of momenta and spins). 
Their expectation values will in general be contaminated by contributions from CP- 
conserving interactions. 

The asymmetries in the t and t polarizations in the production plane, as given in 
eqns. ( |3.6|) and (|3.7|) , translate into T-even observables formed by energies and/or 
momenta of the final states. As an example, we have investigated the expectation 
values of the following two observables (another one was given in |I0|]): 

Ai = E+- E_ (4.1) 

A2 = kt-e^ -kt-r. (4.2) 

Here E±, £^ are the energies and momenta of the leptons in t — > i'^ugb and t i'Uib 
in the laboratory frame and kj(t) is the top (antitop) momentum in this system. 
To measure A2, one has to select events where the t decays leptonically and the 



i hadronically, which in principle allow to reconstruct the i momentum and 
vice versa. We will give explicit expressions for the expectation values of Ai and A2 
below when we discuss contaminations by CP-conserving interactions. In calculating 
these expectation values we have used the narrow width approximation for the top: 
because the top width is much smaller than its mass (in view of the experimental 
upper bound on which is of the order of 200 GeV), the approximation of the on- 
shell production of t and i followed by their weak decays yields a good description 
of the reactions considered here. We also neglected CP violation in the decays of 
t and t (for comments on this, see 0]). For the parton distributions entering the 
calculation of expectation values in pp collisions we have used the parametrization 
of In order to assess the statistical sensitivity of the observables Ai and A2, 



we have computed the signal-to-noise ratios {Ai) / AAi {i = 1,2) for Higgs masses 
100 GeV <m^< 450 GeV both for = 15 TeV (LHC) and = 40 TeV. Here 



AAj = y (Af) — (v4j)2 denotes the width of the distribution of A^. We present our 
results in figs. 7 and 8 for the same parameter set as used in calculating the partonic 
asymmetries, that is, a = —a = l,gvv = 1 (for the definition of gvv see appendix, 
eqn. ( [A.9| )). We integrate here over the whole phase space. Both observables have 



signal-to-noise ratios of order 10 ^. LHC offers larger effects due to the fact that 
AAi 2 is larger for y/s= 40 TeV. 



7 



If both t and t decay leptonically, one can look at the T-odd observable 

T2 = (b - b) ■ {£+ X £-) (4.3) 

where b, b denote the momenta of the b and b jets in the laboratory frame. (This 
observable was also discussed in The expectation value of T2 traces the spin- 

spin correlations of (3^) and ( 3.101 ). In fig. 9 we show the signal-to-noise ratio 



as a function of the Higgs mass for this observable, again with the same choice of 
parameters as for 2. The effect is also of the order of 10~^. 

We will now discuss in some detail possible contaminations of the observables 2 
and T2 due to CP-conserving interactions. Such contaminations arise in particular 
because the pp initial state is not a CP eigenstate. One can give general arguments 
why these contaminations should be small. Most importantly, the dominant subpro- 
cess is gluon fusion which does not induce any CP-conserving contributions to our 
observables (cf. and below). Furthermore, T-odd observables like T2 do not 



receive contributions from CP-invariant interactions at the Born level but only from 
absorptive parts. The main background in this case comes from order and order 
Q^sCtweak absorptive parts in qq ti which generate nonzero functions 63^*^^^ in R^. 
However, numerical simulations show that these contributions are smaller than 10~^, 
i.e. about three orders of magnitude smaller than the signal shown in fig. 9. Poten- 
tially more dangerous are CP-even contributions to Ai and A2, because, as will be 
shown below, they can already be generated by weak interactions at the Born level. 

Integrating over the whole phase space we can actually give explicit analytic for- 
mulae for the expectation values of Ai and A2 in terms of the structure functions. 
This is very illuminating for identifying possible contaminations. We have carried out 
our calculations within the naive parton model (which neglects intrinsic transverse 
momenta of the incoming partons) and restricted ourselves again to gg and qq initial 
states. Furthermore, we have used the narrow width approximation desribed above 
and have taken into account only the SM decays of t and t. Then we find: 



1 1 g{ml,/m^^) 
a 2s 4:71 

I / dxi — £^J-1 / dx2 
^ Jo Xi Jo 

/•I NP{xi) /•! 

+2 / dxi — / dx2 

Jo Xi Jo 



1 , NP{X2) /■! , ^Xi+X2E 



X2 J-1 
NI{X2) rl 
X2 



2y/XlX2 3 ^ 



^g2 



dzf3 



xi + X2 El 



2^X1X2 3 

^^^.(Ei^zA'i + i([(l - z^)mt + z^Ei]W^^^ + E^zW^^^))] }, (4.4) 
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11 g{ml,/mj) 

{A2) -- 



a 2s An 

Xi Jo X2 J-1 L 3 y y 



JO Xi JO X2 J-1 L 3 

- + ^"^twr + ^i^^^D] }■ (4.5) 

4X1X2 o 

Here a is the total cross section for pp — >• ttX, s is the pj9 colhsion energy squared, 
giy) = {1 + 2y + 3y'^)/{2 + Ay), N^, N^(^--^ denote the gluon and quark (antiquark) 
distribution functions of the proton, Ei is the energy of the top quark in the partonic 
CM and 13= {I- ml/ Elf /'^. 

Equations ( [4 .41) and (|4.5| ) exhibit several interesting features: 

- One can see explicitly that gluon fusion generates no CP-even contributions to the 
observables. 

- Quark-antiquark annihilation produces several contaminations: In {Ai) a term 
~ zA'^{z) appears which, after integrating over z, is nonzero only if A'^{z) has a part 
which is odd in z\ that is, if qq ti has a forward-backward asymmetry. Such an 
asymmetry is induced in order a^. (In |l^ this potential source of contaminations 
was discussed.) Possibly more important are the terms 6^™"^ and 6^™"^ which appear 
in both expectation values above because these terms can be generated at the Born 
level via qq ^ Z ^ ti. We calculated their contributions and found that for both 
observables they are suppressed by more than two orders of magnitude in comparison 
to the signals shown in figs. 7 and 8. 

A future multiple TeV and high luminosity collider like the LHC has the potential 
of producing more than 10^ tt pairs. If it were for statistics alone detection of effects of 
a few permil which we found might be feasible. More detailed (Monte Carlo) studies 
including judicious choices of phase space cuts are required in order to explore the 
possiblity of enhancing the signals by some factor. A crucial issue will eventually be 
whether detector effects can be kept at the level of 10~^. 
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5 Conclusions 



In this article we have studied the possibihty of detecting CP violation in top quark 
pair production at future hadron colliders. We have given a general kinematic analysis 
of the underlying dominant partonic subproccsses and identified the relevant CP 
asymmetries at the parton level. We have further computed these asymmetries in 
two-Higgs doublet extensions of the SM where CP violation is generated through 
neutral Higgs boson exchange. Whereas at the parton level these models can induce 
asymmetries of the order of a few percent, realistic observables built up from energies 
and/or momenta of the final states into which the top quarks decay give signal-to- 
noise ratios of up to a few xlO~^. Contaminations by CP-conserving interactions 
were shown to be much smaller than the signals. Since the issue of CP violation is 
of fundamental interest detailed investigations of the experimental feasibility of an 
observation of these effects would certainly be worthwhile. 
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Appendix 



In this appendix we list our analytic results for the structure functions of the tt spin 



density matrix defined in equ. (|2.1| )- (|2.3| ) and decomposed into CP-even and 
CP-odd parts in equ. (p.6|) - (p.8|) and also the corresponding functions in R'^ defined 
in (|2.9|) . All calculations are carried out in the two-Higgs doublet extensions of the 
SM described in section 3. The relevant Feynman diagrams are shown in figures la-h 
for the process gg — > ti and in figures 2a,b for qq — > ti. 

R'^ is obtained very easily: The CP-even part is determined to good approximation 
by the Born diagram fig. 2a, whereas R'cp results from the interference of fig. 2b 
(with couplings aa = —'Jcp) with fig. 2a. The nonzero CP-even structure functions 
of R'^ read: 



CqO 
Cgi 

Cq6 



18 

9tf3\z' - 1) 



9t 
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9 ^{Ei+mt)^ 
-gtP'zE, 



+ 1 



(A.l) 

Here and in the following gs denotes the strong coupling constant, Ei is the energy 



of the top in the CM system of the incoming partons and (3=^1 — m^/Ei. Recall 
that 2 = p ■ k. 

The CP-odd contributions are: 



b^f = ^^^^^^^^^-^^lrnG{s\ (A.2) 



3 



El + nit 



87r2 9 
m^V2Gp^4gtEiPz 
87r2 9 
mlV2GFlcp^9tEif3dz^ - l)f3^Ei 



+ 1 



ImG'(s), (A.3) 



Cql = '''^^Zr^' . ^-Gis), (A.4) 



87r2 9 
Here G^? is Fermi's constant and 



El + mt 



ReG{s). (A.5) 



G{s) - — , (A.6) 
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where 



dH 1 



ZTT^ P — irP +it{l + kiY — mf + ie {I + ki — pi — " f^t + 



(A.7) 



is a standard three-point scalar integral which can be reduced to dilogarithms |^ 
We note here that for the models of section 3 all structure functions are ultraviolet 
finite. In particular, the scalar two-point functions Bq show up in all our results only 
as differences of the form 

B,{qlmlm',)-B,{qlmlml) = 

_|_ x[mi — 1712 " li) + "^2 ~ "^^ 

This completes our results for the matrix W^. 

The computation of is more involved since the contribution of fig. Ih becomes 
resonant if m,^ > 2mt. (Fig. Ih actually represents four amplitudes: two CP- 
conserving ones with couplings and a^, respectively, and two CP-violating ones 
with couplings ad.) The width of ip must therefore be taken into account in the (f 
propagator. We compute by summing the partial widths for ip W~^W~ , ZZ, ti 
in the two-Higgs doublet model which contains (3.1). At the Born level only the 
CP = +1 component of couples to W^W~ and ZZ. The couplings are given by 
the respective SM couplings times the factor 



gvv = {diicos P + d2ism P). (A. 9) 



Explicitly, 



w 



6(m^ — 2mw, 



6(m^ - 2mz) 



glymlV2GFi3w 
glvml^GpPz r 



Pw + 12 



m 



w 



in 



Pi + 12 



rrir. 



6(m^ - 2mt 



?)m^pm'l\/2GFPt 



Stt 



(A.IO) 



1 — Am'^z,t/^%Y^'^- order to incorporate 



Here we have used the notation Pw,z,t 
the resonance region we have determined -Reven from the squared Born amplitudes 
figs, la, lb, the interference of fig. la with the CP-even amplitudes of fig. Ih, and 
the squared amplitudes of fig. Ih. We denote the Born contributions by a lower 
index "Born" and the other two contributions by a lower index "resonance" in the 
following. The results for the nonzero structure functions of -Rfven ^^6- 
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A9 — A9 I A9 
A — ^Born ^ 



resonance 



^Born 



/Iff 
resonance 



192£;f (-1 + /32^2)2V 1 It t 1 t HI 

gi f _ 1 ttvI^Gf mt 

(s-m2)2 + r2m2 I 16 87r2 -1 + /32^2 

[2s(a^/?^ + a2)[ReCo(s, m^^ m^, m^)(s - mj,) 

+ImCo(s, , , ml)V^m^] — Aa^0^{s — m^) 



+ 



m: 



2m2 



32 V 87r2 

+s(a^/?2 + a^a^) |2 - s/?'Co(s, m^, m,^ m^) |2] } 



(A.11) 



resonance 



CgO,Born + CgQ 

,resonancc 

-,K7 + 9/5V) ^(^4_2^2^2^2^4 



192£;f(-l + /32z2)2v--i — I'-t ' 



(s - m2)2 + r2m2 

,2^4 ~2 



{ 



16 87r2 -1 + I3^z^ 



2s(a2/3^ - a^)[ReCo(s, mf, m;)(s - m;) 



+ 



'm3v/2G'i 



^3/ 2~2,o2 ~4\\r^ 2 2 2 



32 V 87r2 

+s(a^/?2 - aV)|2 - s/3^Co(s,mt^mt^mi^)|2]}, 



(A.12) 



g,^/32 (7 + 9/32^2-) (X - ^2 

32(-l + /32^2)2 



(A.13) 



C95 

'^p5,Born 



Cp5,Born ~l~ C(^5, resonance 

gtP' (7 + 9/^2^2) 



-AE^mt^ - 2mt^ - 2(3^Ei^mtZ^ - 2(5^E^^mt^z^ - (5^E^^z^) 
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Cg5, resonance 



1 m.^v^Gi. 



32 



2\ |2 



(A.14) 



96 (El + mt) (-1 + /^^^^ 



,2 • 



(A.15) 



The scalar three point function Co appearing in equations ( [A. Ill ), ( [A. 121 ) and ( [A. 141 ) 
is given by 



dH 1 



1 



ivr^ /2 _ 772^ + (/ — J9i)2 — ml + ie {I — pi — P2) 



ml + ie 



(A.16) 



Numerically, -Reven dominated by the Born contributions. This completes our dis- 
cussion of -Rfven- 

The CP-violating part Rcp results from the interference of the Born diagrams 
with the amplitudes of figs. Ic-lh (with couplings ad = —■Jcp) and the interference 
of the CP-even and -odd amplitudes of fig. Ih. We found that if m^ is of the order 
of 2mt or larger, R^p is dominated in the resonance region by the contributions from 
fig. Ih. Since the complete expressions are rather lengthy, we have split them with 
respect to the contributions from the individual diagrams. For example, means 

the contribution from fig. Ih to the function b'^f and c"^''^^ denotes the part of Cgi 
that is generated by the diagrams of figs. Id and le. 

The function b'j^f gets nonzero contributions only from the box diagrams of fig. 

Ic: 



"si 



''V2Gf7cp 



8n'- 



7 + 9pz)[lmDs{t){l - l3^) - 2lmDu{i)pzEl{l - p^) 
+2(ImDn(t) + lmD2i{i))p^zEl{z^ - ] 
+ (7 - 9(3z)[- lmDs{u){l - f3^) - 2lmDu{u)(3zEf{l - (3^ 
+2{lmDu{u) + lmD2i{u))P^zEl{z^ - 1)] }. 



(A.17) 
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In this expression, 



Ds{i) = Do{i){ml - i) + C^{s,m%mlml) 

Ds{u) = Do{u){m^t -u) + Co{s, mj, ml m^), (A.18) 

(where Co{s,m^,mt,mt) is defined in (|A.7|) ), t = {pi — kiY, u = {P2 — kiY and 

J ZTT'^ r — + 16 [I + ki)^ — mf + ze 
1 1 
{I + ki — pi)^ — ml + ie {I -\- ki — pi — P2Y — mf + ie 

D^it) = Di^{i)ki^ - Di2{i)pi^ - Di^{i)p2^ 

Df,^{i) = D2iii)ki^ki^ + £'22 + D2z{i)p2,,P2v 

-D24ii)ki^Pi^ - D25ii)ki^p2„ + D2&{i)pi^,p2v + D27{i)g^,y. (A.19) 

Dq{u); Dfj_{u); D^^{u) are obtained from ( |A.19|) by interchanging pi and p2- The 
functions Dn, . . . , -D27 can be reduced to expressions which contain only the scalar 
two-, three- and four-point functions Bq, Cq, Dq (see e.g.|^ ). 
The function feS/' reads: 



"g2 

Ac) 

Og2 



ds) 

^92 



'^g2 + '^g2 + "g2 



m: 



'■V2Gf1cp 



{(7 + 9(3 z) \ - ImDsii) ^ ' {E, + mt + (3zE,) 
L -C/i + mt 

+2lmDu{i)mtEi{-Ei + Eiz^ - mtZ^) + A\mD27{i)mt 
+2(Im(L)nt) + \mD2i{i))(3'^El{z^ - l){2mt + Eiz^ - mtz^) 

+ (7 - 9/?^) [ - ImDsiu)^^ — (^1 + mt- (3zEi) 
L Ei + mt 

+2lmDn{u)mtEi{-Ei + Eiz"^ - mtZ^) + AlmD27{u)mt 
+2{lmDn{u) + lmD2i{u))(3^El{z^ - l){2mt + E^z" - mtz")] } 



m: 



'V2Gf1cp ^gtmtf3^zHmG{s) 
87r2 8(-l 
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^92 



m: 



(s - m2 )2 + r2m2 



87r2 4(-l+/32^2) 



(A.20) 



As can be seen seen explicitly from these formulae, all contributions to the functions 
bgf and bg<f result either from absorptive parts of the one loop amplitudes or from 
terms of the form: width times dispersive terms (which is present only in b^^). 
This is in agreement with the general statements made in section 2. One can also 
check the relations following from Bose symmetry as given in table 1. The functions 
Cgl and Cg2 arise from dispersive parts in the one loop amplitude or width terms times 
absorptive parts (which is present only in below). They read: 



Cgl = 



(c) (d,e) (f) 

Si Si Si 



87r2 



96(-l +/52^2) 



{{7 + 9(3z)[ReDs{i){l - - 2RcDn{i)f3zEl{l - (3^) 
-2{ReDii{i) + ReD2i{i))phEl{z^ - 1) 
+(7 - 9(3 z) [ - ReDs(u)(l - (3^) - 2ReDii(u)(3zEl(l - (3^) 
-2(ReDii(u) + ReD2i(u))(3^zEl(z^ - 1)] } 



jd,e) _ mlx^Gplcp 



-9tEi 



87r2 96(-l + /32z2) 



9l3z + 7) [2Co(t, mj, ml ml) {(3^ - 1) 
(3z 



(3z-l 



Cs{t){(3'z'-2(3' + l) 



+(9/3z - 7) 2C^{u, ml, ml ml{p' - 1) 



(3z 



JS) _ 

Cgl — 



Pz + 
m^V2GF-fcp 



-C,(^)(/?V- 2/3^ + 1)]} 



-9t 



87r2 192£;i(-l + /32z2) 



{(7 + 9l3z){Bo{ml mj, ml - Bo{i, mj, ml)^ 



pz-1 



^52 

+{7-9Pz) {Bo{ml m^ ml- Bo{u, ml, ml) 



(3z + 



(A.21) 



where we used the notation 
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and 



C,{i) =Co{i,ml,mlm^,) + 



Boimj, ml, mf) - Bo{i, , mf) 



Cs{u) = CQ{u,m m^,m^) + 



mi — t 

2 _2 , Bajml, m^, mg) - Eo(m, m^, mj) 



ml — u 



(A.22) 



Co(t, ml,ml,ml] 

dH 1 



iTT^ /2 — + ze (/ + /ci)^ — mi + ie{l + ki— piY — mi + ie 



Co{u, mip, mi, mi) is obtained from (|A.23|) by the replacement pi p2. 
Finally, 



(A.23) 



Cg2 
Cg2 



.(d,e) 



(c) (d,e) (f) (g) (h) 
S2 S2 S2 S2 S2 ' 



m: 



7 + 9l3z) 



96(-l + /32z2) 

- ReDs{i) ^ {El + mt~ pzEi) 



El + mt 

+2ReDii{i)mtEi{Ei - Eiz^ + mtz^) + AReD27ii)mt 
+2(ReL>ii(t) + ReD2iii))p^Ei{z^ - l)(2mt + ^1^2- 



mfZ 



+ i7-9Pz) 



ReDJu) 



-{Ei+mt + (3zEi) 



El + mt 

+2ReDii{u)mtEi{Ei - Eiz"^ + mtz^) + 4ReL'27(M)'Tit 
+2(ReDii(M) + ReD2i(ti))/3^^?(-2^ - l)(2mt + Eiz^ - mtz^)] } 



{(7 + 9/5z)[2Co(t,m2,m2,m: 



-Csity 



96(-l+/52^2) 

^ Ei + mt 

Eiz'^ + mt2;2 



(^i+mt-Z^^^i) 



^2(1-/5;^) 



+ (7-9/32;)[2Co(M,m2,m2,m2) 
-mt-Ei^2+mt^ 



- 2m^ - p^Eiz"^ 



El + mt 
2 



{El + mt + pzEi) 



2mi - (3^Ey)] } 
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{(7 + 9(3z)iBoiml , m?) - 5o(t, mj, m^)) 

mi El + rrit — jSzEi 
Ei+mt El{l-pz) 

+(7 - 9/32;)(5o(m2, mj, m^') - 5o(m, , ml)) 

mt Ei + mt + PzEi 
Ei+mt Ef{l + Pz) i 

(g) _ ml^/2Gy^jCi' —3glniij]^z-I{cG{s) 
" 8^i^ 8(-l + p^z^) 

(h) ^ ghcp f _ 1 m\\plGF mtP 

'^^ ~ (s-m2)2 + r2m2 I 16 87r2 -1 + /32^2 

2s(l + /3')[ReCo(s, ml ml ml){s - 
+ImCo(s, , TTij , )r(^m<^] - 4(s - m^) 

+ 1 (^^^^^) ' [rs'PAC.{s, rnl ml m?) ^ 

+2spa^\2 - sP^Go{s,mlmlml\'^]}. (A.24) 
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Table Caption 

Table 1: Transformation properties of the structure functions defined in (2.2)-(2.3). 
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CP 


P 


T 

(ImT = 0) 


CPT 
(ImT = 0) 


"Bose" 


^ l^J 


/I I Z ) 


/I I Z ) 


t\ I Z ) 


A9(y\ 


A9( ■y\ 
/I 1 — Z J 








Ogl\Z) 


Ti^ ( ■y\ 
Ogl[Z) 


( A 
-OgA-Z) 






^'^g2\Z) 


"g2[Z) 


Og2[Z) 


Og2\-Z) 




^gA^) 


Og3[Z) 


( ■y\ 
~0g3[Z) 


-Og3[Z) 


-Og3[-Z) 


CgQ\Z) 


Cgo(^) 


Cgo{z) 


Cgo{z) 




Cgo{—Z) 


C9l(^) 


-Cgl(2:) 


-Cgl{z) 


-Cgliz) 


Cgl(^) 


-Cgl{-Z) 


Cg2(^) 


-Cff2(^) 




-Cgiiz) 


Cff2(^) 


Cg2{-z) 


C53(^) 


-Cff3(^) 




Cgsiz) 


-Cff3(^) 


-Cgsi-z) 


Cg4(2) 


Cg4(^) 




Cg4{z) 


Cg4{z) 


Cg4{-Z) 






C95{z) 


C95{z) 




Cgbi-Z) 






Cgeiz) 


Cgeiz) 


Cp6(^) 


-Cge{-z) 






-Cgviz) 


-Cg7{z) 




Cg7{-z) 


Cg8(^) 


Cg8(^) 


-Cgsiz) 


-Cgsiz) 


-Cp8(^) 


-Cgs{-Z) 



23 



Figure Captions 

Fig. 1: Born level QCD and </? exchange Feynman diagrams which contribute to the 
production density matrix for gg — > tt. Diagrams with crossed gluons are not shown. 

Fig. 2: Born level QCD and relevant ip exchange Feynman diagrams for qq tt. 

Fig. 3: Expectation value (k- (s+ — s_))g as a function of the parton CM energy for 
— 100 GeV (dashed curve) and — 350 GeV (sohd curve). 

Fig. 4: Expectation value (k • (s_|_ x s_))g for the same choice of Higgs masses as in 
Fig. 3. 

Fig. 5: Same as Fig. 3, but for (k • (s+ — s_))g. 
Fig. 6: Same as Fig. 4, but for (k • (s+ x s_))g. 

Fig. 7: Signal-to-noise ratio for the observable Ai (defined in (4.1)) as a function of 
the Higgs boson mass for proton-proton CM energies ^/s — 15 TeV (sohd curve) 
and a/s — 40 TeV (dashed curve). Here = 150 GeV, a — —a — 1, gw — 1- 

Fig. 8: Same as fig. 7, but for the observable A2 (defined in (4.2)). 

Fig. 9: Same as fig. 7, but for the observable T2 (defined in (4.3)). 
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